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Abstract. - We show that the specific heat of the superconductor MgB2 in zero field, for which 
significant non-BCS features have been reported, can be fitted, essentially within experimental 
error, over the entire range of temperature to Tc by a phenomenological two-gap model. The 
resulting gap parameters agree with previous determinations from band-structure calculations, 
and from various spectroscopic experiments. The determination from specific heat, a bulk 
property, shows that the presence of two superconducting gaps in MgB2 is a volume effect. 



The discovery of superconductivity in MgB2 |l| raised the questions of its nature and 
the origin of its relatively high transition temperature T^ ~ 40 K. Specific heat (C) is a 
powerful tool to aid in answering these questions and, more generally, to provide information 
on the thermodynamics of the transition. Several groups have reported such measurements 
on MgB2 [p|-[lO[. It is now established that C significantly deviates from the standard BCS 
behaviour. First, a large excess in C is observed in the vicinity of Tc/4: |g-^. Second, an 
exponential fit of C{T) in the region T ^ Tc indicates a gap ratio 2Ao/kBrc only one- 
quarter to one-third of the isotropic BCS value [^§;^. This excess was interpreted as a 
possible sign of a second superconducting gap, whose existence is predicted by band-structure 
calculations 1 11-13). The specific heat near Tc is puzzling also with the jump AC at Tc 



consistently smaller than the BCS weak-coupling lower bound. In this Letter, we present an 
empirical two-gap model that fits the experimental data over the whole range of temperature 
to Tc. This model resolves the apparent contradiction between different analyses of the specific 
heat, and relevant parameters show good agreement with determinations based on independent 
experiments. 

We focus on two sets of specific-heat data obtained independently in two different lab- 
oratories. Experimental methods and results have been described elsewhere B, S,H, m|. The 
unusual excess specific heat at ~ Tc/4, which denotes the presence of excitations within the 
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main gap, is a consistent feature that is common to different samples and different techniques. 
These measurements also give similar values for the normal-state contribution, with a coeffi- 
cient of the linear term 7n ~ 2.65(15) mJ mol~^ K~^, and satisfy the criterion of the normal- 
and superconducting-state entropy being equal at Tq. However, detailed results, such as the 
height and the width of the jump AC at Tc, are sample-dependent. The sample of Ref. p] 
was a powder of isotopically pure Mg^^B2 embedded in GE7031 varnish, whereas the sample 
of Ref. P] was a sintered commercial powder. A third sample prepared from Mg and B by 
high-pressure techniques gave similar results Q . The electronic part of the specific heat was 
determined by subtraction of the normal-state data, obtained either at fields of 14 or 16 T 
in Ref. , or with a short extrapolation of the 9 T data in Ref. [0 . We refer to the original 
articles for details. 

Although the low-T behaviour of the specific heat data in the earlier studies ^j^ definitely 
pointed to the presence of excitations with a characteristic energy smaller than the BCS gap 
^Bcs = 3.53kBT'c, it was not clear whether this was due to a continuous, but extreme, 
distribution of the gap resulting from anisotropy, or two discrete values of the gap closing at the 
same temperature T^ with possible anomalous temperature dependence at some intermediate 
temperature. Furthermore, it was not clear whether these models could account for the specific 
heat over the whole range of temperature to Tc- We present here a simple empirical model, 
based on the existence of two discrete gaps Ai and A2 at T = 0, both closing at Tc- In order to 
calculate their respective contributions, we first consider the case of a single gap Aq, following 
the method developed by Padamsee et ai, and generally referred to as the a-model flJ. 
The ratio 2Ao/kBT'c (3.53 in the BCS theory) is not fixed, but is considered to be a fitting 
variable. The temperature dependence is taken to be the same as in the BCS theory, i.e. 
A(i) = Ao5{t), where 5{t) is the normalised BCS gap at the reduced temperature t — T/Tc 
as tabulated by Miihlschlegel |13|. The thermodynamic properties, entropy (S) and C, can 
be calculated as appropriate for a system of independent fermion quasiparticles: 



jnTc 7r2 keTc 7o ' Tn^c di 

where / = [exp{f3E) + 1]^^ and f3 — (keT)^^. The energy of the quasiparticles is given by 
E = [e^ + A^(t)]°'^, where e is the energy of the normal electrons relative to the Fermi surface. 
The integration variable is y = e/Ap. 

The fit of experimental data for MgB2 leads to very low values of 2Ao/kBrc for one of 
the gaps, substantially less than 3.53 (see below). The a-model was devised for simulation of 
strong-coupling effects Q, and has usually been applied to strong-coupling superconductors, 
leading to values > 3.53. In that case, the temperature at which the gap closes is lowered 
relative to the normal BCS closing temperature by retardation effects. Since the BCS ratio, 
2Ao/kBTc = 3.53, is the weak-coupling lower limit, smaller values can have no physical mean- 
ing as measures of the strength of the coupling. (However, anisotropy, both as theoretically 
studied [|6| and experimentally observed [l^, does lead to values < 3.5.) In the present case, 
as applied to a two-gap superconductor, a small value of 2Ao/kBTc has no bearing on the 
strength of the coupling , but means only that the temperature at which the small gap closes 
is raised relative to the normal BCS closing temperature by coupling to a larger gap. 

Figure |] shows the calculated C /t^^Tc for 1 < 2Ao/kBTc < 5. We checked the numerical 
results by comparing the data for 2Ao/kBTc — 3.53 with Miihlschlegel's tables |]1^, and 
by verifying that the entropy at Tc is equal to that of the normal state. The curves for 
2Ao/kBrc > 3.5 are similar to those reported in Ref. W%. The unusual shape of the curves 
for low values of 2Ao/kBTc may be understood by considering two characteristic temperatures. 
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Fig. 1 - C/t7nrc vs. t according to the a-model for 2Ao/kBT'c = 1.0, 1.5, 2.0, 2.5, 3.0, 3.5 (thick 
line, BCS), 4.0, 4.5, 5.0. Inset: specific-heat jump at Tc vs. (2Ao/kBrc)^. 



Ta — Ao/(1.76kB), and Tc, which are equal in the BCS limit, but which are independent in 
the present model: 

• For T <C Ta, the thermal energy is too small for many quasiparticles to be excited 
across the gap. Only the tail of the statistical distribution contributes, so that the 
electronic specific heat follows an exponential behaviour approximately, similar to that 
of a semiconductor. 

• Above T « Ta < Tc, the temperature is high enough to excite most of the quasiparticles 
across the gap. The specific heat approaches that of the normal state, although the 
system is still superconducting. 

• At T = Tc, the gap closes. If Tc S> Ta, i.e. if the gap is small compared to the thermal 
energy at Tc, only a small change occurs in the number of excited quasiparticles. The 
BCS ground state is essentially empty. As a consequence, the specific-heat jump is 
small. 

The smaller the gap, the closer the C /t^nTc curve approaches the normal-state line, and 
the smaller the AC at the transition. We verify numerically the relation between the gap 
and the jump, AC = kBiV(0)/(kBTc)2(dAVd/3) ex A^ (inset of Fig. [l]) Hi. This quadratic 
dependence holds only because the variation of the normalised gap with t is common to all 
curves. 

In a two-band, two-gap model, the total specific heat can be considered as the sum of the 
contributions of each band calculated independently according to eq. (|l|) if interband transi- 
tions due to scattering by impurities or phonons can be neglected. Each band is characterised 
by a partial Sommcrfeld constant 7^, with 71 + 72 = 7n- C data are fitted with three free 
parameters, the gap widths Ai and A2, and the relative weights 71/711 = x and 72/711 = 1 — x. 
Figure || shows the data (circles) and the fit (thick line), compared to the BCS specific heat 
(thin line). Insets show the gap functions, and the various contributions to the total electronic 
specific heat. The latter curves show evidence of weak correlation between the fitting param- 
eters; the low-temperature excess is related to A2, whereas the jump at Tc is due essentially 
to the Al component. Numerical results are given in Table ffl. 
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Fig. 2 - BCS normalized specific lieat (thin line), experimental data (o), and two-gap fits (thick 
lines), versus the reduced temperature t. (a) data from Ref. M; (b) data from Ref. |2|. Insets: gaps 
2Ai/kBrc and 2A2/kBrc versus t (dotted lines), and partial specific heat of both bands (full lines) 
Parameters obtained from the fits are given in Table 
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Fig. 3 - Semi-logarithmic plot of the electronic specific heat versus 1/T. Dashed line: asymptotic 
curve, eq. (^j with 7 = Tn; thick line: eq. (g) with 7 = 0.47n (see text); thin line: standard BCS 
curve, also shown in Fig. H; (o), data from Ref. M. 

Fig. 4 - Superfiuid fraction versus reduced temperature. Thin line: contribution of Ai; dotted 
line: contribution of A2; thick line: full two-gap fit; (o): data obtained from measurements of the 
penetration depth presented in Ref. psi . Fitted parameters are given in Table IE 
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Table I - Gap ratios 2Ai/kBT'c, 2A2/kBTc and weights x as determined by the two-gap model (lines 
1-4) o.nd by different techniques (lines 5-10). 

Ref. Technique 2Ai/kBrc 2A2/kBrc x : [1 - x) 



: 1 specific heat 


4.4 


1.2 


55% 


45% 


\ 1 specific heat 


3.8 


1.3 


50% 


50% 


: 1 specific heat 


3.9 


1.3 


50% 


50% 


25 penetration depth 


4.6 


1.6 


60% 


40% 


19 Raman 


3.7 


1.6 






20 photoemission 


3.6 


1.1 






21 tunneling 


4.5 


1.9 






22 point-contact spectroscopy 


4.1 


1.7 






23 point-contact spectroscopy 


4.2 


1.0 






13 band structure 


4.0 


1.3 


53% 


47% 



In spite of its limitations, this empirical model fits the measured specific heat well over 
the whole range of T to Tc. The sample dependence of the results is reasonably low, and may 
reflect metallurgical differences. The larger value of Ai for the sample of Fig. |p (isotopically 
pure Mg^^B2 powder [0) reflects a sharper jump and a steeper slope just below T^ compared 
to the sample of Fig. |b (MgBz sinter (|]). On average, 2Ai/kBrc -- 4.0 and 2A2/kBrc - 1.2, 
with approximately equal weights. 

Moreover, the fitted parameters are qualitatively and quantitatively comparable with inde- 
pendent determinations from other sources. They are consistent with band-structure calcula- 
tions |13[ and spectroscopic measurements [19^E3[ , which report the presence of two gaps, the 
smaller gap having approximately one-third the BCS value and the larger gap being slightly 
greater than the BCS value (Table M. We emphasise that C, a thermodynamic property, 
probes the whole volume, whereas spectroscopic measurements are more sensitive to surface 
conditions. 

The relative weights (1:1, i.e. x ^ 0.5) are consistent with the calculations of Ref. |l^. Liu 
et al. attribute the larger gap Ai to particular 2D sheets of the Fermi surface, whereas the 
smaller gap A2 is associated with 3D sheets. Using partial densities of states and de Haas-van 
Alphen mass renormalizations, the weight of the smaller gap is evaluated as a; ~ 0.47, and 
1 — cc ~ 0.53 for the larger one. The agreement with the two-gap model fits is remarkable. 

The present two-gap model reconciles the apparently conflicting results of Ref. [g and |^,^ . 
By fitting their specific heat data close to Tc, Kremer et al. [g| concluded that their data was 
consistent with a medium- to strong-coupling 2Ao/kBT'c ~ 4.2. However, the fitted value of 7n 
at Tc was 1.1 mJ mol^^K"^, less than half of 7n measured in the normal state. Alternatively, 
Yang et al. Q and Bouquet et al. Q fitted the exponential decrease of the low-T data and 
concluded that 2Ao/kBrc ^ 0.9. However, the fitted value of 7n at low T was too small also, 
0.7 mJ mol~^K~^ in Ref. g. In the framework of the two-gap model, the main contribution 
just below Tc is that of the larger gap Ai, with a break in the slope characteristic of medium- 
to strong-coupling, and an amplitude of AC determined by 71 = xj^ ~ 7n/2 (insets of Fig. 0), 
in qualitative agreement with Kremer's analysis. The main contribution at T <IC Tc is that of 
the smaller gap A2, with the exponential decrease determined by A2, and the amplitude by 
72 = (1 — x)7n ~ 7n/2 (insets of Fig. 0), again in qualitative agreement with the analysis of 
Ref. ^ . The latter data are presented below in a slightly different approach. Rather than the 
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usual empirical interpolation C oc exp(— 1.44Tc/T), we use the low-T asymptotic formula [M: 

In Fig. 0, we plot data in the form ln(C/7nr) versus 1/T, together with the limit given by 
eq. (||). With 2Ao/kBTc = 0.9 and 7 — 711, eq. (g) overestimates the data, although the slope 
determined by Aq is correct. With 7 = 0.47n, the fit is good in the domain where eq. (||) 
holds. 

The same two-gap model can be applied to the superfluid density p, which is given, for 
a single gap, by p = 1 — 2Ao/kBT/p°° /(I — f)dy. The penetration depth A oc p~^/^ is 
given in Ref. |23] and is plotted in Fig. 0, together with a two-gap fit (thick line) and its 
components (full and dotted lines). These data are not strictly bulk measurements, but 
probe the sample to a typical depth of A = 1800 A p|. Nevertheless, A is large compared to 
the typical sampling depth of many spectroscopic experiments, which is on the scale of the 



coherence length ^ = 50 A 
other determinations (Table 



The fitted parameters 2l^i/\i^Tc and x are consistent with 

!)• 

The empirical a-model allows a quantitative comparison to be made between different 
experiments and theory within a general framework. The results are numerically consistent, 
and confirm the coexistence of two gaps for the bulk sample. This situation holds the promise 
of interesting single-crystals properties. Our two-gap model is phcnomcnological since we 
postulate the existence of the gaps, without specifying their origin. Any theoretical approach 
leading to a similar average electronic density of states would be compatible with the present 
results, so that specific-heat measurements alone cannot settle in favour of any particular 
microscopic model ||ll|,[l3|. 

Some limitations exist. First, the a-model assumes a BCS-like T-dependence of the 
gap. However, if the variation of the smaller gap is reasonably smooth, the results should 
not depend critically on its exact shape, since the main effect on the specific heat occurs 
below Ta where A2(T) is expected to be essentially constant. Self-consistent calculations of 
A(r) might lead to corrections, and more elaborate simulations are currently under way ||2q| . 
Second, we calculate each contribution of the gaps independently and assume that they are 
additive. Some coupling is present, but within the present model, its sole effect amounts to 
bringing the natural closing temperature of the smaller gap, i.e. « 10 K, up to « 40 K. 

Our two-gap model describes only the zero-field specific heat. As data in iJ > suggest a 
different field dependence for each gap, a theory of the mixed-state specific heat for a two-gap 
superconductor would be most useful in extracting quantitative information from C(T,H). 
Indeed, the field dependence of the electronic contribution at low temperature is unusual. 
The coefficient of the linear term in the mixed-state, "f{H), dramatically increases in small 
fields |0-^, in a quasi-logarithmic way g, and saturates for fields much below Hc2, in fact 
near i/c2/2. Moreover, the characteristic dip in C/T for < T < 10 K associated with the 
gap A2 in one of the bands vanishes by « 0.5 T. Qualitatively, these results would seem to 
indicate that a small field is able to quench the smaller gap, in agreement with spectroscopic 
measurements p2|. The saturation of 7(i?) much below iJc2 suggests that a major part of the 
electrons are either normal or in a gapless state, possibly by virtue of inter-band scattering 
in the presence of a normal sheet on the Fermi surface. Forthcoming models may have to 
embody the different dimensionality of the gaps. Interesting developments are expected for 
the physics of the vortex state for superconductors with such an unusual fc-dependent gap. 
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